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Abstract 

A time-dependent wave-packet approach is presented for the quantum dynamics study of the 
AB-fCDE reaction system for zero total angular momentum. A seven-degree-of-freedom calculation 
is employed to study the chemical reaction of H 2 +C 2 H H 4- C 2 H 2 by treating C 2 H as a linear 
molecule. Initial state selected reaction probabilities are presented for various initial ro- vibrational 
states. This study shows that vibrational excitation of H 2 enhances the reaction probability, 
whereas the excitation of C 2 H has only a small effect on the reactivity. An integral cross section 
is also reported for the initial ground states of H 2 and C 2 H. The theoretical and experimental 
results agree with each other very well when the calculated seven dimensional results are adjusted 
to account for the lower transition state barrier heights found in recent ab initio calculations. 


Electronic address: dywang@nas . nasa . gov 



INTRODUCTION 


2 


• Developments of quantum dynamics scattering calculations 

— First quantum dynamics calculation on H+H 2 system 0 

— First 3D quantum dynamics calculation on H+H 2 in 1976 6 

— First 6D quantum dynamics calculation on H 2 +OH in 1994 c 

• During the past two decades, quantum scattering theories well developed 
for 3 and 4 atom systems 

• Now dynamics calculaions compete with experiment in their accuracy 

— H+D 2 d and F+H 2 e give the same results as experiments 

— Some properties of H+H 2 0, D 2 CF 


“ F.Z. London, Elektrochem. 35, 532 (1929); H. Eyring, and M.Z. Polani, Phys. Chem. (Munich), B12, 279 
(1931). 

b G.C. Schatz and A. Kuppermann, J. Chem. Phys. 65, 4642 (1976). 

c D.H. Zhang and J.Z.H. Zhang, J. Chem. Phys. 101, 1146 (1994); D. Neuhauser, J. Chem. Phys. 100, 
9272 (1994). 

d M.P. de Mranda, D.C. Clary, J.F. Castillo and D.E. Manolopoulos, J. Chem. Phys. 108, 3142 (1998). 
e D.E. Manolopoulos, J. Chem. Soc., Faraday Trans. 93, 673 (1997). 
f D.H. Zhang, M.A. Collins and S.-Y. Lee, Science 290, 961 (2000). 


INTRODUTION 


3 


• More than 4 atom systems - quantum dynamics computations still a 
formidable task. 

— almost all the calculations are reduced dimensional 

— a full dimensional quantum dynamics study for the H+CH^ 

— no full dimensional quantum dynamics been done on 5 atom systems 

• Why study H 2 +C 2 H? 

— C 2 H is one of the most abundant interstellar polyatomic molecules de- 
tected^ 

— found in Titan and Jupiter’s atmospheres* 

— dominant chain carrier to form acetylene 

— hydrocarbon reactions provide disequilibrium species 

— only 4D calculations available for this reaction J ' k 

— PES of Wang and Bowman^ 


9 F. Huarte-Larranaga and U. Manthe, J. Chem. Phys. 113, 5115 (2000); J. Chem. Phys. 116, 2863 
( 2002 ). 

h K.D. Tucker, M.L. Kuther and P. Thaddeus, Astrophys. J. L115, 193 (1974). 

7 D.F. Strobel, Planet. Space Sci. 30, 839 (1982); G.R. Gladstone, M. Allen and Y.L. Yung, 119, 1 (1996). 
9 D. Wang and J.W. Bowman, J. Chem. Phys. 101, 8646 (1994). 

*H. Szichman, M. Gilibert, M. Gonzalez, X. Gimenez, and A. Navarro, J. Chem. Phys. 113, 176 (2000); 
ibid , J. Chem. Phys. 114, 9882 (2001). 



THEORETICAL METHODS 


• Jacobi coordinates for the reaction AB+CDE 


Full-dimensional Hamiltonian 

H = ~ + (J Z0 1 + hn(ri) + K{r2) + K[r3) + 
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THEORETICAL METHODS 


ID vibrational reference Hamiltonians 

hri( r i) = ~ 2 ^ ,^ r 2 VlD( r i)i * — 1 , 2 , 3 
Total angular momentum function of J23=j2 + J*3 in BF frame 


y;gf 23 (r 2 , f 3 ) = Y D & L(*> ^ &)■ 


2 j 2 + 1 


< j3«T-i 2 0|jr'2377Z > Y h m(6 3 ) 


Overall internal angular momentum rotational function of the AB+CDE 
system, j 123 = ji + j 23 




Y < ilT — ^23^23^23 Iil23-fir > 0l) 

^23 

x Yi 2 f 23 (r 2 , f 3 ) 


TD wavefunction expanded in terms of the BF rovibrational eigenfunctions 


iT f J Me 
y -vjK 


nviv 2 vzhhhK ( 5 ) 

* Y™AnK{R,h,h,h) 


THEORETICAL METHODS 


Relationship between the SF rotational basis function ^ ^2, ^ 3 ] 

and BF rotational function Y/^j^kA n, r 2 , h) 


Y Zh* K foh,h,h) 


2L + 1 


2 J + 1 


< ji23^ L0\J K > 


r[,r'2,r' 3 ) 


The SF rotational function 2// 1 ^ 3 j 123 i(-R',^I, f' 2 , r' 3 ) 


Vhhu iaiC#. *i> *£) = E < iusmjLM - m,| JM > 


my mi m 2 


X < jimij^rrij - mi\j m m,j > 
x < j2m 2 j 3 m z \j 2z m j - m 1 > 

x 2/i 3 (^ 3 ) 2/72 (^2 ) 2/71 (n ) 2/x (-R) 
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Parity operator e on the BF rotational basis function gives 


iY hhumK(R, h ,f 2 , f 3 ) = e £ 


2L + 1 


2J + 1 


< J123K L 0 \J K > 


xAmL&SlAA) 

= J2(-l) jm+L+J < jm _ 7O,0| J - > 


x {-l) h+i2+h+L y™u mL (R' > f[,r' 2 , r' 3 ) 

= (-^ 1+j2+h+jm+JY nn 3 h, >- k(R, h,h, r 3 


• Parity adapted BF total angular momentum eigenfunction 

Y /JriwK(R,h,r2,r 3 ) 


Y hhzimi<{R, h, h, h) = y 2(1 + fl ’ f2 ’ ^ 

+ f x , f 2 , f 3 




SEVEN-DIMENSIONAL TREATMENT 


• 7D Hamiltonian 

n 2 d 2 , (J-ji 2 ) 2 


H = - 


+ 


+ 


2fidR 2 2 fiR 2 


+ hri( r l) + ^r 2 ( r 2) + ^r 3 ( r 3) + 


Jl' 


2/xirf 


32 


(10) 


2 wf 


+ V 7 d (R, r u r 2 , r 3 , 6>i , <9 2 , 0i) 


^2? ^ 3 ) becomes two angular momentum coupling case 


Y &!k w =J^rh^) [Y ^ K{9u 02 > ^ 


(ii) 


+ 2 _*(0i, * 2 , <h)] 


• 7D initial wave-packet 


= 0) = $*o(#)^i 0 ( r i)^W r2 ’ r ^ Y hXh2 0 Ko ( 12 ) 


• ^^ 23 0 i 2 0 ( r 2 , 7 * 3 ) is the eigenfunction of the two dimensional vibrational Hamil- 
tonian 


H = h r2 (r 2 ) + ^ 3 ( 7 * 3 ) + V 2D (R = 00 , r u r 2 , r 3 , 0i, 0 2 , 0 3 = 0, </>i, </> 2 ) (13) 




SEVEN-DIMENSIONAL TREATMENT 
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• Split-operator method to propagate the wave packet 

'ip JMe ( R, n, r 2 , r 3 , t + A)= e~ lH ° A/2 e~ lUA e~ lH ° A/2 'ip JMe (R, r 1 , r 2 , r 3 , t) 

( 14 ) 

• Reference Hamiltonian Hq defined as 

Ho = + M r 0 + + ^ 3 (^ 3 ) (15) 

• Effective potential operator 


U = 


(J-J12) 2 , jl 


2fiR 2 


+ 


2 /xi r 2 


+ 


J2 

2/^2 r\ 


+ V 7D {R , n, r 2 , r 3 , 0i, 0 2 , 0i) (16) 
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RESULTS AND DISCUSSION 


• Initial state reaction probabilities for ji=0) + 0211(^2=^3=0, j2=0) 

— y H + C 2 H 2 as a function of kinetic energy for J=0 




Reaction Probability 


RESULTS AND DISCUSSION 
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• Initial state reaction probabilities for H2(yi=0, ji) + C2H(i/2=^3=0, j 2 =0) 
— > H 4- C 2 H 2 as a function of kinetic energy for J=0 
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RESULTS AND DISCUSSION 

• Initial state reaction probabilities for H2 (Vi= 0, ji= 0 ) + C2H(z/2=^3=0, j2) 
— > H + C2H2 as a function of kinetic energy for J =0 
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• Comparison between present 7D resulsts k?D, rescaled 7D results kjjj x 
e 2.08 jkT an( ^ x e 2.24/kT us j n g barrier height differences between Wang 

and Bowman’s PES and ab initio results^, and various experimental results 
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1 X. Zhang Y. Ding, Z. Li, X. Huang and C. Sun, J. Phys. Chem. 104 , 8375 (2000); Y. Kurosaki and T. 
Takayanagi, J. Chem. Phys. 113 , 4060 (2000). 
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• Conclusions 

— vibrational excitation of H 2 molecule enhances the reaction probabili- 
ties;while the vibrational excitation of C 2 H doesn’t change the reactivity 
much 

— ground rotational state of H 2 gives the largest probability;^ = 1 initial 
rotational state of C 2 H gives the largest probability 

— Both the 7D quantum calculation results and the ab initio studies 
indicate the current barrier height of Wang and Bowman’s PES is too 
large 

• Future work 

— Extend the current 7D calculation to full-9D calculation 

— Extend the AB+CDE calculation to A+BCDE calculation 




